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NASA’s first Air Traffic Management (ATM) Technology Demonstration (ATD-1) was 
created to facilitate the transition of mature ATM technologies from the laboratory to 
operational use. The technologies selected for demonstration are the Traffic Management 
Advisor with Terminal Metering (TMA-TM), which provides precise time-based schedul- 
ing in the terminal airspace; Controller Managed Spacing (CMS), which provides terminal 
controllers with decision support tools enabling precise schedule conformance; and Interval 
Management (IM), which consists of flight deck automation that enables aircraft to achieve 
or maintain a precise spacing interval behind a target aircraft. As the percentage of IM 
equipped aircraft increases, controllers may provide IM clearances to sequences, or strings, 
of IM-equipped aircraft. It is important for these strings to maintain stable performance. 
This paper describes an analytic analysis of the string stability of the latest version of 
NASA’s IM algorithm and a fast-time simulation designed to characterize the string per- 
formance of the IM algorithm. The analytic analysis showed that the spacing algorithm 
has stable poles, indicating that a spacing error perturbation will be reduced as a function 
of string position. The fast-time simulation investigated IM operations at two airports 
using constraints associated with the midterm airspace, including limited information of 
the target aircraft’s intended speed profile and limited information of the wind forecast on 
the target aircraft’s route. The results of the fast-time simulation demonstrated that the 
performance of the spacing algorithm is acceptable for strings of moderate length; however, 
there is some degradation in IM performance as a function of string position. 
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I. Introduction 


O VER the next twenty years, the Federal Aviation Administration (FA A) is predicting a substantial increase 
in the number of revenue passenger miles flown. 1 This increase could strain the current air traffic system 
if improvements are not made. To help increase the efficiency of arrivals, new procedures, such as Optimized 
Profile Descents (OPDs), have been implemented at a number of sites. However, these procedures are often 
not used during periods of peak traffic due to the lack of precision scheduling and spacing tools, 2 resulting 
in sub-optimal arrival operations. The National Aeronautics and Space Administration (NASA) has been 
developing technologies that are expected to increase the use of efficient arrival procedures during periods 
of peak traffic. 

The goal of NASA’s first Air Traffic Management (ATM) Technology Demonstration (ATD-1) is to 
accelerate the implementation of mature NASA technologies, enabling aircraft to use speed control to remain 
close to their OPDs during periods of peak traffic. 3 Three NASA technologies were selected to be integrated 
together to achieve this goal: the Traffic Management Advisor with Terminal Metering (TMA-TM), 4,5 which 
provides precise time-based scheduling in the terminal airspace; Controller Managed Spacing (CMS), 6,7 
which provides terminal controllers with decision support tools that enable precise schedule conformance; 
and Interval Management (IM), which consists of flight deck automation that enables aircraft to achieve 
or maintain a precise spacing interval behind the preceding aircraft. 8,9 These three technologies extend 
precision scheduling into the terminal area and provide pilots and controllers with control mechanisms that 
improve schedule conformance and in-trail spacing precision. 

The IM algorithm that NASA has been developing is a trajectory-based IM algorithm called the Airborne 
Spacing for Terminal Arrival Routes (ASTAR) algorithm. 10, 11 The ASTAR algorithm uses the intended 
trajectories of both the IM and target aircraft to calculate their Estimated Times of Arrival (ETAs). The 
ETAs are used to calculate a spacing error and the speed control required to null the spacing error by an 
achieve- by point. Previous versions of the ASTAR algorithm have shown good performance in several batch 
simulations and human-in-the-loop simulations. 12 18 

The ATD-1 flight demonstration will occur in a midterm airspace environment where future commu- 
nication infrastructure, such as controller-pilot data link communication, is not available. To compensate 
for the lack of accurate target aircraft trajectory intent information, NASA’s ASTAR spacing algorithm 
uses constraints from published Standard Terminal Arrival Routes (STARs) to compute the ETAs of both 
the IM and target aircraft. Discrepancies between the published speed profile and delayed speed profiles 
computed by TMA-TM cause the speeds flown by the target aircraft to deviate from the published speeds 
expected by the spacing algorithm. Additionally, the IM aircraft is not expected to have access to the target 
aircraft’s wind forecast for the ATD-1 demonstration, increasing the uncertainty of the target aircraft’s ETA 
calculations when the target aircraft is on a different route. 

The latest version of the ASTAR algorithm, ASTAR- 12, includes modifications to the speed control 
algorithm that improve compatibility with TMA-TM and CMS when a delayed target aircraft deviates from 
its expected groundspeed. A previous batch simulation examined pairs of aircraft, consisting of a single IM 
aircraft following a single target aircraft, conducting IM operations using the ASTAR- 12 algorithm. 19 The 
simulation found that the ASTAR- 12 algorithm’s delivery accuracy and precision were generally acceptable; 
however, the batch simulation did not examine the performance of the ASTAR- 12 algorithm when there 
is a string of several IM aircraft following each other. This paper describes both the fundamental string 
stability of the ASTAR- 12 spacing algorithm and a fast-time simulation designed to emulate strings of 
aircraft performing IM operations in the midterm national airspace system. 

The first section of this paper provides an overview of previous IM research. The following section 
describes the fundamental string stability behavior of the ASTAR- 12 algorithm using an analytic approach. 
The final two sections describe a fast-time simulation that investigated the performance of strings of IM 
aircraft using the ASTAR-12 spacing algorithm, assuming midterm airspace constraints. 

II. Background 

A. Previous IM Algorithm Research 

IM research has been conducted since the 1970’s. The early IM work focused on the development of IM 
algorithms, such as time history algorithms and constant time delay algorithms, that used state information 
from the target aircraft to determine the spacing error. 20 Within these algorithms, the spacing error is 
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defined as the difference between a spacing goal assigned by air traffic control and the time elapsed between 
the time when the target aircraft passed the IM aircraft’s current position and the time when the IM aircraft 
arrives at that position. Due to the definition of spacing error, these algorithms are limited to operations 
where the IM and target aircraft are in-trail. 

More recently, trajectory-based IM algorithms have been developed. The NAS A- developed Airborne 
Spacing for Terminal Arrival Routes (ASTAR) is a trajectory-based IM algorithm that was designed to be 
used in a NextGen environment where detailed information of the target aircraft’s intended trajectory is 
available to the IM aircraft. 11, 12 The ASTAR algorithm uses predicted trajectories of the IM aircraft and 
target aircraft to calculate a time-to-go for each aircraft to a common achieve- by point. The spacing error 
is then calculated by comparing the difference between the time-to-go of the IM and target aircraft with a 
spacing goal assigned by air traffic control. Several batch simulations and human-in-the-loop simulations 
have demonstrated that the ASTAR algorithm is able to precisely deliver aircraft to the achieve-by point 
and that the speed produced by the algorithm are generally acceptable to pilots. 8,9, 17,21 Each of these 
human-in-the-loop simulations assumed an advanced airspace environment with controller-pilot datalink 
communications used to transmit IM clearances to the flight deck, and that the IM aircraft will have access 
to detailed information of the target aircraft’s intended trajectory. 

When the ATD-1 project began, the focus of IM research at NASA switched from a future environment 
that included controller-pilot datalink communications to the use of IM in the midterm airspace environment. 
Since controller-pilot data link communication is not expected to be available in the midterm national 
airspace system, IM clearances are provided using voice communications and the intended trajectories of the 
IM and target aircraft are assumed to be published Standard Terminal Arrival Routes (STARs). 3,15,16,22 In 
preparation for the ATD-1 flight demonstration, simulations were also conducted to examine the integration 
of IM with TMA-TM and CMS (also known as Terminal Sequencing and Spacing). 

The previous version of the ASTAR algorithm, ASTAR-11, did not perform well when used with TMA- 
TM and the CMS tools. TMA-TM uses the predicted trajectory of each aircraft along their projected OPDs 
to compute their ETAs to a series of scheduling waypoints. If there is a conflict at one of the scheduling 
waypoints, TMA-TM often delays aircraft. With the advent of flex scheduling, aircraft can also be advanced 
in certain circumstances. Since ASTAR used the published STARs as the estimate of the target aircraft’s 
intended trajectory, the speeds expected by ASTAR and TMA-TM do not always match. The ASTAR-11 
algorithm was not designed to compensate for large speed differences between the target aircraft’s actual 
speed and the published speeds that were used by ASTAR-11 to predict the target aircraft’s ETA. The result 
was that the ASTAR-11 algorithm exhibited a large steady state error and undesirable closure rates with 
the target aircraft when the target aircraft was absorbing delay. 19 

NASA’s ASTAR algorithm was updated to mitigate the previously described problems and improve com- 
patibility with TMA-TM and CMS. 19 The main modification was a groundspeed term that was added to 
the ASTAR algorithm to compensate for discrepancies between the target aircraft’s actual speeds and pub- 
lished speeds. The groundspeed term essentially enables the IM aircraft to match the target aircraft’s speed 
deviation and then correct for the spacing error using the proportional control term. The groundspeed term 
also prevents steady-state errors from occurring when the the target aircraft is not flying its expected speed, 
reducing undesirable closure rates between the IM and target aircraft. Both human-in-the-loop simulations 
and batch simulations investigated the performance of the ASTAR- 12 algorithm and its acceptability to 
pilots and controllers. However, these studies did not focus on the string behavior of the ASTAR-12 spacing 
algorithm. 

B. String Stability Research 

When there are a number of IM-equipped aircraft arriving into an airport, it is possible for several IM 
aircraft to follow each other, creating a string of IM aircraft. A spacing algorithm is considered to be string 
stable if the spacing error does not grow as a function of string position, ensuring that small perturbations 
early in the string do not develop into unacceptable behavior later in the string. 

Figure 1 shows an example of a string of IM aircraft. The aircraft are depicted as being in-trail; however, 
when a trajectory-based algorithm is used, aircraft within a string can either be on in-trail or on different 
routes. Throughout the analytical analysis described in section III, increasing values of i refer to moving 
backwards in the string and decreasing values refer to moving forward in the string. 

There is a substantial amount of research that has investigated the string stability of various IM algo- 
rithms. Weitz and Hurtado examined the string stability of several IM algorithms using both an analytic 
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Figure 1. Example of a string of IM aircraft 

approach and a fast-time simulation. 23 One of the algorithms they examined was a time-to-go algorithm 
that is similar to a previous version of the ASTAR spacing algorithm (ASTAR-11). They found that the 
time-to-go algorithm was weakly string stable. Krishnamurthy et al. conducted a fast-time simulation that 
examined strings consisting of 100 IM aircraft. 24 The results of the fast-time simulation did not reveal any 
string instabilities. 

The string stability of various spacing algorithms have also been validated in previous human-in-the-loop 
simulations. In 2005, Barmore et al. described a human-in-the-loop simulation that included strings of nine 
IM aircraft. The results of the simulation did not reveal any string instabilities. 13 Baxley et al. conducted a 
human-in-the loop simulation that examined a previous version of NASA’s ASTAR spacing algorithm. 18,21 
During the simulation, strings of 23 IM aircraft, including six pilot flown aircraft, were examined for string 
instabilities; no string instabilities were found. 

C. ASTAR-12 Spacing Algorithm Description 

ASTAR- 12 is a trajectory-based spacing algorithm designed to support IM operations both when the IM 
and target aircraft are in-trail and when they are on merging routes. The spacing error is calculated using 
the time-to-go of the IM and target aircraft along their predicted 4D-trajectories. In the future, it will be 
possible for either the target aircraft’s trajectory or estimated time of arrival at the achieve- by point to be 
communicated directly to the IM aircraft. However, this information is not available in the midterm airspace 
environment. To compensate for the lack of 4D-trajectory information, published Standard Terminal Arrival 
Routes (STARs) are used to predict both the IM and target aircraft’s 4D-trajectories. The time-to-go for 
each aircraft is simply the difference between their ETA at the achieve-by point and the current time. The 
spacing error is then defined as the difference between IM aircraft’s time-to-go ( TTGi ) to the achieve-by 
point and the target aircraft’s time-to-go to the achieve-by point (TTGi- 1) minus the spacing goal assigned 
by air traffic control (A). 


ei(t) = TTGi(t ) - TTGi-i(t) - A (1) 

The amount of speed control required to null the spacing error is computed using a proportional con- 
troller with a groundspeed term added to compensate for differences between the target aircraft’s predicted 
groundspeed and actual groundspeed. The speed control is added to the IM aircraft’s nominal speed to 
generate an IM commanded speed 

Vcmdi (t) = Vnorrii (t) T k p ei(t) T kGs{vd eri _ 1 (t) Vnorrii-i W) (2) 

Here v nomi is the nominal 4D-trajectory airspeed of the IM aircraft; k p is the proportional gain, e*(t) is 
the spacing error; kcs is the gain for the groundspeed term; Vden-% is the airspeed derived from the target 
aircraft’s groundspeed, wind forecast, and altitude; and v norni _ 1 is the nominal 4D-trajectory airspeed of the 
target aircraft. 

The IM avionics used in the ATD-1 flight demonstration will likely be a retrofit implementation that is 
not connected with the aircraft’s autoflight system, requiring the flight crew to monitor for IM speed changes 
and enter them into the aircraft’s mode control panel speed window to close the control loop. To reduce the 
number of speed changes that the flight crew is required to respond to, speed commands are discretized into 
either five or ten knot increments prior to being displayed to the flight crew. 

Since the ASTAR-12 algorithm is designed to support an achieve-by IM operation, the ideal behavior is to 
smoothly null the spacing error between the point where the IM operation begins and the achieve-by point. 
Gain scheduling is implemented to ensure that the appropriate amount of speed control is used throughout 
the arrival. The proportional gain ranges from a value of 0.375 when the IM aircraft is far from the achieve- 
by point to 1.5 when the IM aircraft is within 10 nmi of the achieve-by point. The groundspeed gain is 
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set to one until the IM aircraft is 40 nmi from the achieve- by point. Between 40 nmi and 20 nmi from the 
achieve- by point, the groundspeed gain is linearly decreased to zero, and remains at zero until the achieve- by 
point. Additionally, the groundspeed gain is set to zero whenever the target aircraft’s groundspeed is faster 
than expected. This is done to reduce the number of instances where ASTAR-12 commands a speed increase 
during an arrival and to increase conformance with the controller tools, which primarily resolve conflicts at 
meter points during high density operations by delaying aircraft. 

Additional filtering is applied to the target aircraft’s groundspeed differential to prevent undesirable speed 
changes from occurring. The target aircraft’s groundspeed is filtered using a first order low pass filter to 
remove high frequency variability. The time constant of the filter changes from 60 seconds when the IM 
aircraft is more than 35 nmi from the achieve-by point to zero when the IM aircraft is at the achieve-by 
point. 


III. Analytic Analysis of String Stability 

Long strings of IM aircraft could exhibit unstable string behavior if the control algorithm is not designed 
properly. This section describes an analytic analysis of the string stability of the ASTAR-12 algorithm. The 
first step of the analysis is to derive a transfer function for the error propagation through a string of IM 
aircraft. After the transfer function is derived, its poles can be calculated and used to determine if the error 
propagation transfer function is asymptotically stable. As is often the case with analytical analyses, several 
simplifying assumptions are made within this analysis: 

• The time delay associated with various filters in ASTAR is not modeled; 

• The time it takes the flight crew to recognize a speed command and enter it into their aircraft is not 
modeled; 

• The conversion from airspeed to groundspeed is not modeled; 

• The nominal speed of the IM and target aircraft are assumed to be both constant and equal; and 

• Simplified aircraft dynamics are used. 

The string stability of the ASTAR-12 control algorithm can be computed similarly to the analysis of 
the time-to-go algorithm described by Weitz and Hurtado. 23 A simple one dimensional point mass is used 
to model the longitudinal aircraft dynamics. The aircraft’s acceleration is assumed to be proportional to 
the difference between the IM commanded speed (v cm d) and the IM aircraft’s current speed (t^), with a 
proportional gain equal to k v . The variable x refers to the distance traveled by the aircraft and the subscript 
i refers to the i th aircraft in the string. 


Xi(t) = Vi(t) 


( 3 ) 


0 — k v (v C md{t) ( 4 ) 

In a time-to-go algorithm, such as ASTAR-12, the spacing error is defined as the difference between 
the time-to-go of the IM aircraft and the time-to-go of the target aircraft minus the assigned spacing goal 
(equation 1). Assuming that the nominal profile speed is constant and ignoring the conversion from airspeed 
to groundspeed, the time-to-go is simply the along-path distance to the achieve-by point ( xabp ~ x(t)) 
divided by the nominal profile speed ( v nom ). Thus, the spacing error and its first and second order time 
derivatives can be described using equations 6, 7, and 8. 


ei(t) 


ei(t ) 

e’i(t) 


xabp Xj(t ) xabp - ^ 

Xnom X norn 

Xj-i(t) - Xj(t) _ 

Xnom 

- Vj(t) 

Xnom 

- Vj(t) 

Xnom 


( 5 ) 

(6) 

( 7 ) 

(8) 
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The ASTAR-12 speed control law described in equation 2 is combined with the aircraft dynamics described 
in equation 4. 


k v (v norn 4“ k p Ci{t^ T ^GS 1 (^) ^ nom ) (9) 

= k v (v norn + fe p ei_i(t) + k GS {vi- 2 (t ) - i; nom ) - v*_i(t)) (10) 

Next, equations 9 and 10 are used in equation 8, and the resulting equation is rearranged and simplified 
into a differential equation describing the error propagation through the string of IM aircraft. 


ei(t) = 


kv(v n0 m T kp6{— l(t) T kQs(vi—2{t) Vnom) 1(^)) 

Vnom 

ky {Vnom T k p 6i{t) T kQs(v<i—\(t) Vnom) ^z(^)) 


( 11 ) 


kyVnom k v V nom ^ k v k p ei-i(t) k v k p ei(t) k v k G sv nora ^ k v k G sv nora 

Vnom Vnom V nom V nom V nom V nom 

+ k„fc S ( _ kv ( 


( 12 ) 


k v k p ei-i(t) k v k p €{(t^ 


+ k v k GS ei-i(t) - k v ei(t ) 


(13) 


The first and second terms in equation 12 cancel each other, the fifth and sixth terms cancel each other, 
and the seventh and the last term are equal to e^_i(t) and e^(t), respectively. After simplification, equation 
13 is reduced to: 


&i(t) + k v ei(t ) + 


= k v k G sei-i(t) + 

'Vnom 


k v k p 


Vnom 




(14) 


A transfer function for spacing error is calculated by taking the Laplace transformation of equation 14 
and rearranging. 


Ei - 1 k v k G ss + 

L - 1 - u nom 

Ei s 2 +k v s+E *he. 

u nom 


(15) 


The result is a second order transfer function with two poles and one zero. Within control theory, a 
system is considered asymptotically stable if the real portion of all the poles are negative. When the poles 
are negative real numbers, the system response will be a damped exponential, and if the poles are complex 
conjugates with a negative real portion, the response will be a damped oscillation. Equation 15 can be 
rewritten into an augmented version of the standard canonical form that includes a zero in the numerator. 


Ei -1 = (P + I)^n 

Ei s 2 + 2C w„s + u 2 


(16) 


Here, uj n is the natural frequency, £ is the damping ratio, and — z is the zero (i.e., the value of s where 
the numerator of equation 16 equals zero). In a standard second order system without any zeros, damping 
ratios greater than one will result in exponential decay, a damping ratio between zero and one will result in 
a damped oscillation, and a damping ratio of zero will result in a marginally stable osculation that does not 
grow or decay. A negative zero results in increased rise time (i.e., faster control response), and increased 
overshoot than the standard system. 


h h 

rv v rv p 


(17) 
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The poles of the transfer function in equation 16 are determined by finding the values of s where the 
denominator is equal to zero. 


poles = v n (-( ± \j C 2 -* 1) (20) 

It can be shown that the real portion of the poles of the transfer function are less than zero if k v , k pi 
and v nom are positive, indicating that the ASTAR-12 algorithm is string stable. Additional intuition on the 
effect of the gains on the string stability can be gained by examining equation 18. Increasing values of k p and 
decreasing values of v nom will result in a lower damping ratio. During an arrival, the values of k p increase as 
the aircraft approaches the achieve- by point and the nominal profile speed ( v norn ) decreases, indicating that 
the ASTAR-12 algorithm will have the smallest damping ratio when the aircraft is close to the achieve- by 
point. 

The zeros of the transfer function in equation 16 are determined by finding the values of s where the 
numerator is equal to zero. As the zero approaches — oo, the effect of the zero diminishes and the system 
behavior approaches the standard second order canonical form. Since the value of v norn is significantly larger 
than the value of fc p , the effect of the zero on the string stability of ASTAR-12 is small. 

zero = — - — — (21) 

kGS^nom 

Figure 2 shows the step response for a variety of damping ratios. Assuming that typical values of k v 
range from 0.05 to 0.2, values of k p range from 0.375 to 1.5, and values of v norn range from 100 knots to 
600 knots (the lowest and highest groundspeed expected by a turbojet), the range of damping ratios of the 
ASTAR-12 algorithm can be determined. Using this information, the damping ratio can range from a value 
of 0.9 to a value of 8.9. 



Figure 2. The step response of the error propagation function for several 
damping ratios (k v =0.1, fees — 1> and v norn / k p is varied) 


IV. Fast-Time Experiment Design 

A fast-time simulation was conducted to evaluate the performance of strings of aircraft conducting IM 
operations using the ASTAR-12 spacing algorithm. The simulation examined 6,720 strings of IM aircraft 
arriving into either Phoenix Sky Harbor airport (KPHX) or Denver International airport (KDEN). Each 
IM aircraft conducted IM operations from a point approximately 150 nmi from the runway threshold to 
the final approach fix, which was used as the achieve-by point in this simulation. The experiment variables 
were selected to simulate realistic variability within the midterm national airspace system. The speed 
profiles of the first target aircraft (the leader of each string of IM aircraft) at KPHX were directly derived 
from aircraft speed profiles observed in an ATD-1 human-in-the-loop simulation and wind fields observed at 
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KPHX and KDEN were used. Other aircraft parameters were randomly selected to more accurately simulate 
the variability expected in future IM operations in the midterm National Airspace System. 

Since the objective of this simulation was to characterize the string performance of the ASTAR-12 spacing 
algorithm, the string position was the independent variable of interest. Three additional scenario- level 
variables were included in the simulation to add variability. These variables were: the speed profile flown by 
the first target aircraft in each string, the truth and forecast wind fields, and the destination airport (figure 
3). 

There were a total of 14 target aircraft speed profiles simulated for each airport. The goal when creating 
the target aircraft speed profiles was to approximate the speeds that a delayed target aircraft would fly 
if controlled by an air traffic controller using TMA-TM and the CMS. The target aircraft speed profiles 
investigated at KPHX were based off of speed profiles observed in a previous ATD-1 human-in-the-loop 
simulation (figure 4). To ensure that a wide range of target aircraft speed profiles were included in this 
simulation, the target aircraft speed profiles observed in the human-in-the-loop simulation were clustered 
to determine unique patterns. Twelve speed profiles were selected from the clustered dataset and combined 
with two nominal target aircraft speed profiles. Since human-in-the-loop simulation data was not available 
for KDEN, 12 off-nominal target aircraft speed profiles were created with the goal of approximating the 
variability of the KPHX target aircraft speed profiles. The 12 off-nominal target aircraft speed profiles were 
combined with two nominal speed profiles. 

A total of sixteen wind cases, eight for each airport, were selected from a set of winds observed at KPHX 
and KDEN. Each wind case consisted of a truth wind and a forecast wind that were temporally offset in 
time. The ETA differences between more than a years worth of observed winds and a no wind condition 
were computed for the routes simulated in this experiment, and used to select a set of truth wind fields with 
a range of ETA differences. Wind fields within six hours of each selected truth wind were used as potential 
forecast wind candidates. The Root-Mean-Square (RMS) of the wind difference between each chosen truth 
wind and each forecast wind candidate were examined, and the forecast wind was selected to have a RMS 
difference of 10 to 12 knots in the center airspace and 7 to 9 knots in the terminal airspace. To simulate the 
decent wind forecast that is normally available on-board the aircraft, the forecast wind field was discretely 
sampled using a vertical column at the runway threshold at the altitudes of 30,000 feet, 24,000 feet, 18,000 
feet, 10,000 feet, and at the surface. This method of sampling the forecast winds is consistent with the 
capability that will likely be available in retrofit IM avionics. The discrete descent wind forecast was used by 
the ASTAR algorithm to predict the time-to-go of both the IM and target aircraft. The ASTAR algorithm 
also includes a heuristic that blends the IM aircraft’s sensed winds with the wind forecast; a functionality 
that was used in this simulation. 

In addition to the previously discussed variables, several aircraft parameters were randomly varied. To 
approximate conditions observed in ATD-1 human-in-the-loop simulations, the initial delay of each aircraft 
was selected from a normal distribution with a mean of 30 seconds and a standard deviation of 15 seconds. 
The aircraft dynamics model was randomly selected from a set of seven Base of Aircraft Data (BAD A) 25 


Phoenix Sky Harbor Airport (KPHX) 


Denver International Airport (KDEN) 




Figure 3. Experiment Design 
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(a) East flow EAGUL arrival at KPHX 




(c) North flow BOSSS arrival at KDEN 



Figure 4. The target aircraft speed profile options 


aircraft models using a uniform distribution. The aircraft models used in this simulation include the B737, 
B752, B773, A320, A319, A306, and CRJ1 models. The route of each aircraft in each string was randomly 
selected from a set of two routes at KPHX and two routes at KDEN using a uniform distribution, creating 
different merging and in-trail conditions. The route options that were simulated at KPHX were the EAGUL 
and MAIER arrivals in the east flow configuration to runway 08, and the route options that were simulated 
at KDEN were the BOSSS and ANCHR arrivals in the north flow configuration to runway 35L. Lastly, the 
weight of each aircraft was randomly selected from a set of weights that ranged from the aircraft’s reference 
weight to the aircraft’s reference weight plus half of the maximum payload weight. A total of thirty replicates 
for each combination of independent variables were used to sample the random variables (i.e., 30 replicates 
of each cell in figure 3). 

In addition to the randomized variables, a 10 second time delay was implemented to model the time 
it takes for pilots to recognize a new IM commanded speed and manually enter it into their aircraft. The 
10 second time delay is similar to the flight crew time delay determined from previous human-in-the-loop 
simulations. 21,26 

Aircraft state data from both the target aircraft and IM aircraft were collected and used to compute 
several metrics that describe the performance of the ASTAR-12 algorithm. These metrics include the delivery 
accuracy and precision at the achieve- by point, the schedule drift, the number of commanded speed changes, 
the Root-Mean-Square (RMS) of the spacing error, and the amount of speed control used. 
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V. Fast-Time Simulation Results 


A. Delivery Accuracy and Precision 

The delivery accuracy and precision are measures of how accurately and precisely the ASTAR-12 algo- 
rithm achieved the desired spacing goal. The delivery accuracy was measured as the difference between the 
assigned spacing goal and the elapsed time between when the target and IM aircraft crossed the final ap- 
proach fix. Previous experiments have shown the ability of IM algorithms to achieve precise spacing intervals 
with both a mean delivery accuracy and standard deviation less than five seconds. 13,16,19,21,26 Overall, the 
mean delivery accuracy in this simulation was 2.4 seconds and standard deviation was 4.6 seconds. A total 
of 4.2% of the IM aircraft had delivery accuracies worse than 10 seconds. 

Parsing the data by airport revealed that there was a significant difference in both the mean and standard 
deviation of the delivery accuracy. The mean delivery accuracy at KPHX was 2.1 seconds with a standard 
deviation of 3.0 seconds, and the mean delivery accuracy at KDEN was 2.7 seconds with a standard deviation 
of 5.7 seconds. Additionally, only 0.8% of the aircraft at KPHX had delivery accuracies worse than 10 seconds, 
whereas 7.5% of aircraft at KDEN had delivery accuracies worse than 10 seconds. The large discrepancy 
between the delivery accuracies at KPHX and KDEN is likely caused by differences in the wind conditions 
or the target aircraft speed profiles used in this simulation. 

The delivery accuracy was also analyzed to determine if it changed as a function of string position. 
An Analysis of Variance (ANOVA) and Tukey multi-comparison test revealed that the delivery accuracy 
increased as a function of string position at both KPHX (p <0.001) and KDEN (p <0.001). However, the 
delivery accuracy results did not reveal any unstable string behavior, and the 12th aircraft in the string had 
an average delivery accuracy of 3.3 seconds and a standard deviation of 3.8 seconds, which is within the 
acceptable range. While the delivery accuracy degraded as a function of string position, figure 5 suggests 
that the increase reached an asymptote at the seventh or eighth aircraft in the string. 

Another observation was that the delivery accuracy had a late bias after the fourth aircraft in the string. 
In the next section, it will be shown that the bias in the delivery accuracy resulted in an increase in the 
average schedule drift as a function of string position. While there is no conclusive information on the cause 
of the delivery accuracy bias, the most probably cause is the lack of symmetry in the groundspeed feedback 
term. The groundspeed term in the ASTAR-12 algorithm is disabled when the target aircraft’s groundspeed 
is faster than predicted to prevent undesirable speed increases from occurring. Based on the examination of 
case studies, random variations throughout the string can cause an IM aircraft to be commanded a speed 
faster then the nominal profile close to the achieve- by point. When this occurs, the groundspeed term is 
turned off and the ASTAR-12 algorithm reverts to proportional control. The combination of a target aircraft 
speed deviation and proportional control causes the ASTAR-12 algorithm to exhibit a steady-state error; 
potentially resulting in the delivery accuracy bias. 

B. Schedule Drift 

The schedule deviation, or schedule drift, was also examined as a function of string position. Large 
schedule drifts could cause integration problems with non-IM aircraft that are controlled to achieve a sched- 
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Figure 7. A boxplot of the schedule drift at the final 
approach fix as a function of string position 
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Figure 8. The modeled schedule drift at the final ap- 
proach fix, adjusted for the observed schedule deviation 
bias, compared to the actual schedule drift (error bars 
represent ± one standard deviation) 


tiled time of arrival. Since there were no scheduling tools modeled in this simulation, the schedule drift was 
computed as the sum of spacing errors of all of the aircraft in each string. 

Figure 7 shows that both the mean and standard deviation of the schedule drift increased as a function 
of string position. The mean schedule drift increased from -0.2 seconds for the first aircraft in the string 
to 28.9 seconds for the last aircraft in the string. Likewise, the standard deviation of the schedule drift 
increased from 6.7 seconds for the first aircraft in the string to 28 seconds for the last aircraft in the string. 
The change in the mean schedule drift is caused by the late bias in delivery accuracy, whereas, the increase 
in standard deviation is an expected result based on a model described by Levitt. 27 

Levitt et al. derived an equation that can be used to estimate the schedule drift of a string of IM aircraft, 
assuming normally distributed spacing errors with zero mean. 27 The equation uses the standard deviation 
of the delivery accuracy ( a spc ) and the standard deviation of the first target aircraft’s schedule deviation 
(&sch) to model the standard deviation of the schedule deviation ( cr dev ) for the k th aircraft in the string. 

a dev = + k a spc (22) 

Using the standard deviation of the IM aircraft (4.6 sec) and assuming that the standard deviation of 
the schedule deviation of the first target aircraft in the string is zero, the theoretical standard deviation of 
the schedule drift for the 12th aircraft in the string is 16 seconds. The observed standard deviation of the 
schedule drift for the 12th aircraft in the string was 30 seconds. The theoretical calculation assumes a mean 
spacing error of zero and thus a mean schedule drift of zero; however, it is easy to adjust the model to include 
a delivery accuracy bias. The blue line in figure 8 shows db one standard deviation of the modeled results 
after they were adjusted for the observed delivery accuracy bias (i.e., 2.4 seconds was added to the mean 
for each string position). The larger standard deviation of the observed data is caused by the fact that the 
delivery accuracy was not normally distributed. Instead, several outliers conditions contributed to longer 
tails on the delivery accuracy distribution. 

C. Number of Speed Changes 

The number of speed changes commanded by the ASTAR-12 spacing algorithm were evaluated to ensure 
that they did not increase as a function of string position, which could indicate string stability problems. 
Figure 9 shows that the number of IM speed commands increased slightly as a function of string position; 
the first IM aircraft in the string had an average of 10.6 IM speed commands and the 12 th IM aircraft in the 
string had an average of 12.9 IM speed commands, including published speed changes. This corresponds to 
an average speed change rate of 0.4 speed changes per minute with a standard deviation of 0.1 speed changes 
per minute throughout the arrival for the first aircraft in the string, and 0.5 speed changes per minute with 
a standard deviation of 0.1 speed change per minute for the 12 th aircraft in the string. 

A second metric that was examined was the number and magnitude of speed reversals commanded by 
the ASTAR-12 algorithm. The number and magnitude of speed reversals were computed by counting the 
number of consecutive speed increases. For instance, if there were three 10 knot speed increases in a row 
followed by a speed decrease, one 30 knot speed reversal was counted. Figure 10 shows the number of speed 
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Figure 9. The number of speed changes commanded by 
the ASTAR-12 spacing algorithm as a function of string 
position (error bars represent zb one standard deviation) 


Figure 10. The number and magnitude of speed rever- 
sals commanded by ASTAR-12 as a function of string 
position (error bars represent zb one standard deviation) 


reversals as a function of string position, binned by the speed reversal magnitude. There was an average of 
2.3 speed reversals per IM operation; however, approximately half of the speed reversals were 10 knots or 
less. 

Figure 10 shows that the total number of speed reversals increased slightly as a function of string position. 
More interestingly, the average number of 10 knot speed reversals decreased as a function of string positions 
and the number of higher magnitude speed reversals (greater than 20 knots) increased. 

D. Root-Mean- Square of Spacing Error 

The main indication of whether an IM algorithm is string stable is whether the spacing error increases 
as a function of string position. The RMS of the spacing error over the entire IM operation was examined 
to validate the results of the analytic model presented in section III. Since each string position was provided 
with the same distribution of initial spacing errors (figure 13), an increase in the RMS of the spacing error 
as a function of string position could indicate unstable string performance. Figure 11 shows that the RMS of 
the spacing error did not increase as a function of string position, validating the analytic analysis conducted 
in section III. 

E. Root-Mean-Square of Speed Control 

The RMS of the speed control throughout an arrival was examined to determine if the total amount 
of speed control required to achieve the assigned spacing goal increased as a function of string position. 
Since each string position was provided with the same distribution of initial spacing errors (figure 13), an 
increase in the total amount of speed control as a function of string position could indicate undesirable string 
performance. An Analysis of Variance (ANOVA) and Tukey multi-comparison test revealed that there was a 
substantial increase in the total amount of speed control as a function of string position at KPHX (p< 0.001 ) 
and at KDEN ( p<0.001 ). Figure 12 shows that the average RMS of speed control for the first IM aircraft 
in the string was 11.4 knots and increased to 17.9 knots for the 12th aircraft in the string. The increase in 
the amount of speed control required is likely the driver behind the small increase in the number of speed 
changes and the increase in the number and magnitude of speed reversals. This data suggests that the speed 
commands of aircraft at the end of a very long string of IM aircraft will be less desirable than the commands 
for the first aircraft in the string. 


VI. Conclusions 

Both an analytic analysis and a fast-time simulation were conducted to characterize the string per- 
formance of the ASTAR-12 spacing algorithm. The analytic analysis included the derivation of an error 
propagation transfer function. Both the poles of the transfer function and the damping ratio indicated that 
the ASTAR-12 spacing algorithm is string stable. The fast-time simulations examined 6,720 strings of IM 
aircraft arriving into either KPHX or KDEN. The independent variables that were examined were the des- 
tination airport, string position, wind field, and the speed profile of the first target aircraft in each string. 
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RMS of Spacing Error (sec) 



String Position 

Figure 11. The Root-Mean-Square (RMS) of the spacing 
error as a function of string position (error bars repre- 
sent ± one standard deviation) 


RMS of Speed Control (knots) 



Figure 12. The Root-Mean-Square (RMS) of the speed 
control used to achieve the spacing goal as a function 
of string position (error bars represent ± one standard 
deviation) 


Initial Spacing Error (sec) 



String Position 

Figure 13. The initial spacing error for each aircraft in 
the string was sampled from the same distribution (error 
bars represent ± one standard deviation) 


Additional aircraft parameters were randomly varied to better simulate the variability expected in future 
IM operations. The data was examined to characterize the string performance of the ASTAR-12 spacing 
algorithm. The results did reveal degradation of ASTAR-12 performance as a function of string position; 
however the performance is still in the acceptable range for strings of moderate length. 
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